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^ ■ Abstract. 

Our Galaxy is a barred spiral. Recent work based on the COBE NIR 
i-" ' data implies a small bulge-bar and a disk with a short scale-length. The 

[t . corotation radius of the bar is in the range 3-4.5 kpc. The stellar density 

distribution beyond the end of the bar appears to be perturbed strongly 
by the Galaxy's spiral arms. 

Gas flow calculations in corresponding potentials provide a qualita- 
tive explanation of many features observed in HI and CO Zti-diagrams. 
K^ \ These include the 3-kpc-arm and the apparent four-armed spiral struc- 

^^ ■ ture between corotation and the solar radius. 

^sq \ The mass of NIR-luminous matter is constrained by the terminal 

^SJ ■ velocity curve, the Oort limit, and the bulge microlensing observations, 

^D \ and this implies that the Milky Way has a near-maximum disk and a 

2^ ' dark halo with a large core radius of ~ 15 kpc. 

However, we are still some way from a detailed quantitative model 



l> 



a\ 



i-Q ' for the large-scale dynamics of the Galaxy. I summarize a number of 



o 
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uncertainties as well as how future work might resolve them. 



1. Introduction: The Barred Milky Way 



\^ . It is now well-established that the Milky Way is barred. This is an important 

^ \ development: some long-standing questions about Galactic dynamics are begin- 

ning to be answered in the new framework, as discussed in later sections, and, 
not least, it is changing the way in which we have to think about the Galaxy's 
evolutionary history. 

The strongest evidence for a rotating bar in the inner Galaxy comes from the 
NIR light distribution, source count observations, the atomic and molecular gas 
morphology and kinematics, and the large optical depth to microlensing. This 
section gives a brief summary of this evidence; other recent reviews concerning 
observations and dynamical implications of the Galactic bar are in Gerhard 
(1996), Kuijken (1996), and Morris & Serabyn (1996). 

Observations of cold gas in the inner Galaxy reveal large non-circular mo- 
tions (e.g.. Burton & Liszt 1978 [HI], Dame et al. 1987 [^^CO], Bally et al. 1988 
[^^CO]). Some of the more prominent features include the 3-kpc-arm, the 135- 
km/s-arm, the molecular parallelogram or 180-pc-ring, and the high central 
peak in the terminal velocity curve at / ~ ±2°. Many papers in the past have 
suggested that the observed kinematics are best explained by gas motions in a 



barred potential (Peters 1975, Cohen & Few 1976, Liszt h Burton 1980, Ger- 
hard & Vietri 1986, Mulder k Liem 1986, Binney et al. 1991, Wada et al. 1994). 
In §5 I discuss new hydro dynamical studies of the gas flow in the barred Milky 
Way. 

The COBE-DIRBE photometry clearly show (Weiland et al. 1994, Freuden- 
reich 1998, see Fig. 1) that the Galactic bulge is both brighter and more extended 
in latitude at given positive longitude than at the same negative longitude, ex- 
cept for a region close to the Galactic Centre where the first effect is reversed. 
These signatures are just as expected for a triaxial bulge with its long axis in the 
first quadrant (Blitz & Spergel 1991). Detailed modelling results are discussed 
in §2. 

The flux distribution for IRAS bulge sources is systematically brighter on 
the I > side (Nakada et al. 1991), arguing for a non-axisymmetric bulge. A 
similar signature is also found in the much larger OGLE clump giant sample 
(Stanek et al. 1997). The small intrinsic luminosity spread (~ 0.2 — 0.3 mag) 
makes these stars good distance indicators and allows to constrain the param- 
eters of the bar. Nikolaev & Weinberg (1997) reanalyzed the IRAS variable 
population, for which distance information is available from the known range of 
AGB star luminosities, and again find a barred source distribution. 

The inferred optical depth in microlensing experiments towards the bulge 
(Udalski et al. 1994, Alcock et al. 1997) exceeds the values predicted by ax- 
isymmetric mass models by a factor of 2 — 3 (Kiraga & Paczynski 1994, Evans 
1994). This has widely been taken as further evidence for the Galactic bar. As 
discussed in §3, this is probably correct, but the precise interpretation of the 
measured event rates needs to be clarified and the predictions from models in- 
corporating the bar need to be improved before quantitative conclusions can be 
drawn. 

This article is not a comprehensive review of the 'Dynamics of the Galaxy', 
but is restricted to a number of topics related directly or indirectly to the ex- 
istence of the bar. In the following I will discuss recent work on the structure 
of the bulge and disk (§2), bulge microlensing (§3), the stellar dynamics of the 
bulge §4), the gas flow in the Galactic disk interior to the Sun (§5), and the 
mass of the bulge and disk including consequences for the distribution of dark 
matter in the Galaxy (§6), and end with some concluding remarks in §7. 

2. Photometric structure of the Galactic bulge and disk 

The currently best models for the distribution of old stars in the inner Galaxy 
are based on the COBE/DIRBE NIR data. The COBE data have complete 
sky coverage, and provide broad-band emission maps from J,K,L,M in the NIR 
to 100 /i and 240// in the FIR. However, their analysis is complicated by the 
relatively low spatial resolution, the residual effects of dust absorption, and the 
fact that they contain no distance information; the best results will ultimately 
be obtained by combining them with other complementary data sets. Because 
extinction is important towards the Galactic nuclear bulge even at 2//m, the 
first task is to correct (or 'clean') the DIRBE data for the effects of extinction. 
Arendt et al. (1994) did this by assuming that the dust lies in a foreground 
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Figure 1. Left: Asymmetry map for the cleaned COBE L-band data of 
Spergel et al. (1996). Contours are spaced by 0.05 mag; dotted contours 
indicate larger flux on the I < 0-side. Right: Asymmetry map of a model 
obtained by Lucy- Richardson deprojection, assuming a bar angle oi (j) = 25° . 
From Bissantz et al. (1997). 



screen, while Spergel, Malhotra &: Blitz (1996) and Freudenreich (1998) used 
fully three-dimensional models of the dust distribution. 

The cleaned DIRBE data confirm that the Galactic bulge is both brighter 
and more extended in latitude (b) at positive longitudes than at corresponding 
negative longitudes (/), except for a region close to the Galactic Centre where 
the bulge is brighter at negative / (Dwek et al. 1995, Bissantz et al. 1997). These 
signatures are just as expected for a triaxial bulge with its long axis in the first 
quadrant (Blitz & Spergel 1991). They can be quantified by an 'asymmetry 
map' - flux at (/ > 0,6) divided by flux at {—l,b) in logarithmic units - cf. 
Fig. 1 (from Bissantz et al. 1997). The region of negative asymmetry at small |^| 
argues for a bar rather than a lopsided distribution as the cause of the positive 
asymmetry at larger |/|; see also Sevenster (1997). Most of the asymmetry signal 
is in the range 0.1 — 0.3 mag. 

Binney & Gerhard (1996) developed a Lucy-Richardson deprojection algo- 
rithm to interpret these data, based on the assumption that the Galaxy's NIR 
emissivity distribution is eight -fold (triaxially) symmetric . They showed that 
in this case the ambiguity inherent in the deprojection of 2d data is reduced 
to essentially that involved in choosing the orientation of the bulge's symmetry 
planes. This algorithm was used by Binney, Gerhard & Spergel (1997, hereafter 
BGS97) to fit non-parametric models for the 3d emissivity j{r) to the cleaned 
L-band data of Spergel et al. (1996). When the orientation of the symmetry 
planes is fixed, the recovered emissivity j(r) appears to be essentially unique, 
but physical models for the COBE bar can be found for a range of bar orienta- 
tions: 15° ^ (/> ^ 35°, where (p measures the angle in the Galactic plane between 
the bar's major axis at I > and the Sun-centre line. Zhao (1999) has given 
an illustration of the non-uniqueness of the bar orientation in terms of the even 
part of the bulge density distribution. 

Fig. 2 shows the deprojected luminosity distribution resulting for = 20°. 
This shows an elongated bulge with axis ratios 10:6:4 and semi-major axis ~ 
2kpc, surrounded by an elliptical disk that extends to ~ 3.5 kpc on the major 
axis and ~ 2 kpc on the minor axis. There is a maximum in the projected NIR 
emissivity ~ 3 kpc down the minor axis in the Galactic plane, which is caused by 
emission near the plane and is probably due to incorrectly deprojected strong 





Figure 2. Symmetrized luminosity model for the inner Galaxy obtained by 
Lucy-Richardson deprojection of the cleaned COBE L-band data for (/> = 20°. 
Left: Density projected along the z-axis, contours spaced by 0.1 dex. Right: 
Isodensity surfaces in the zx and zy planes, contours spaced by 0.2 dex. Axis 
lengths in kpc. From Binney, Gerhard & Spergel (1997). 



spiral arms (see below). This feature corresponds to the ring-like structure 
discussed by Kent, Dame &: Fazio (1991). An exponential fit at larger radii 
gives a radial disk scale length of Rd = 2.5 kpc. Fig. 1 shows the asymmetry 
map of a similar model to compare with the data. 

Dwek et al. (1995) and Freudenreich (1998) have constructed parametric 
models to match their dust-cleaned DIRBE data. In Freudenreich's (1998) 
recent work the region near the Galactic plane is excluded in the fit (owing to 
uncertainties in the dust extinction) and the disk is modelled with a central 
hole. His analysis based on several NIR wavebands gives a bulge with axial 
ratios 10:4:3, length 3.1 kpc, at a preferred orientation of around (j) = 14°. The 
disk scale-length outside the central hole is Rd = 2.6 kpc. 

The differences between these models perhaps give a fair measure of the 
remaining uncertainties. It is likely that the dominant source of error is still 
buried in the details of the dust absorption model, but also systematic errors in 
constructing 3d models from 2d data have to be explored further. One major 
difference is in the structure of the inner disk and, hence, the amount of lumi- 
nosity in the inner disk compared to the bulge. To resolve this clearly requires 
fitting the data at small I and h and hence a reliable dust model. Absorption 
will be a lesser problem at higher latitudes where the bulge structure is mainly 
determined, but the observed asymmetry signal of a few tenths of a magnitude 
is not greatly larger than the estimated errors in the extinction (~ 0.07mag, 
Spergel et al. 1996). Thus it appears worthwhile to reanalyze the DIRBE data 
based on more sophisticated dust extinction models (e.g., Sodroski et al. 1997). 

A different and very promising approach is based on analysing large point- 
source catalogues. Stanek et al. (1997) fitted bulge-bar models to the apparent 
magnitude distributions of clump giant stars in 12 OGLE fields. The small 
intrinsic luminosity spread {^ 0.2 — 0.3 mag) makes these stars good distance 
indicators. Stanek et al. 's best bulge model fitted to the reddening-corrected 
data is one with an exponential density distribution, axial ratios 10:4:3, and a 
bar angle of (/) = 20 — 30°. These parameters stem from a parametric fit of the 



bulge only and some contamination from the foreground disk seen in the NIR 
is likely. With potentially many more fields to come from the OGLE II and 
MACHO projects this method has future potential. 

Nikolaev & Weinberg (1997) reanalyzed the IRAS variable population, for 
which distance information is available from the known range of AGB star lu- 
minosities. This sample has a large spatial coverage near the plane; taking into 
account selection effects and incompleteness, Nikolaev & Weinberg (1997) find 
a barred distribution (dominated by the disk population) with length ~ 3.5 kpc, 
axial ratio in the plane 10:4, and orientation ~ 21°. 

NIR starcounts have confirmed the longitudinal asymmetry due to the bar 
(Unavane & Gilmore 1998), at a level more in accord with the BGS97 model. 
Lopez-Corredoira et al. (1997) have modelled the Two Micron Survey Starcounts 
in three strips across the bulge. Their best model is a nearly end-on {(j) = 
12° ±6°) triaxial bar with axis ratios 10:5-6:3. Earlier starcount models (Robin 
et al. 1992, Ortiz & Lepine 1993) also favour a short disk scale {Rd — 2.5 kpc), 
but the uncertainties due to inbuilt model assumptions are not easily quantified. 
Earlier claims that the Milky Way disk could have a central hole have not been 
confirmed (Kiraga, Paczynski & Stanek 1997). 

In summary, this combined work has confirmed and elucidated the barred 
structure of the Galactic bulge and inner disk. While the length of the bar 
is near 3-3.5 kpc and the bar angle (p with respect to the Sun-Center line is 
probably around = 20 — 25°, there is still some debate about the axis ratios 
(within a range 10: [4-6]: [3-4]) and the detailed density structure. A further 
important conclusion with considerable consequences for Galactic mass models 
(see §8) is that the radial scale-length of the Galactic disk in the NIR is short, 
Rd — 2.5 kpc. Further progress will in my view come from (i) improved models 
for the dust extinction, (ii) combining the area coverage of the COBE data with 
the distance information in the clump giant and other stellar samples, and (iii) a 
better understanding of possible population differences between bulge and disk. 

3. Microlensing towards the Galactic bulge 

Microlensing observations provide new constraints on the structure of the Galac- 
tic bulge and disk. The most robust observable is the total optical depth av- 
eraged over the observed fields. Defining r_6 = r/10~®, the measured values 
are: for OGLE main-sequence stars, r_6 = 3.3 it 2.4 (2cj) from 9 events (Udalski 
et al. 1994); for MACHO main-sequence stars, r_6 = 1.9 ± 0.8 {2a) from 41 
events (Alcock et al. 1997), when a correction for disk contribution is removed; 
and r_6 = ?>.S^'^2a (^^) ^^r MACHO's 13 clump giant stars, averaged over ~ 10 
square degrees centered at (/,6) = (2.55°, —3.64°) (Alcock et al. 1997). 

The interpretation of these numbers is complicated. It has recently become 
clear that blending or amplification bias (Alard 1997; blended sources below the 
survey magnitude limit become bright enough to be included when lensed) is im- 
portant for the bulge microlensing events. About half of all OGLE events show 
strong centroid shifts (Goldberg & Wozniak 1998), hence are affected by blend- 
ing. The required correction to the optical depth has yet to be determined. The 
effect is likely to be similarly important for the MACHO main sequence stars, 
whereas the bright clump giant stars should be unaffected. A proper analysis 



of the larger MACHO 4 yr dataset is eagerly awaited. Another complication 
are the steep gradients in the bulge expected from the NIR models; spatially 
resolved data would of course be ideal for structural information. 

It has long been known that axisymmetric models predict r_e ~ 1 — 1.2, 
insufficient to explain the quoted optical depths (Kiraga & Paczynski 1994, 
Evans 1994). Models with a nearly end-on bar as described in §2 enhance r 
because of the longer average LOS from lens to source. The maximum effect 
occurs for (p ~ arctan(6/a) when rbar/''"axi — {sin2(f))~'^ ~ 2 for cf) = 15° (Zhao 
& Mao 1996). In addition, r increases with the mass and the length of the 
bar/bulge. 

Optical depths in Baade's window in some recent bar models for the Milky 
Way are r_6 ^ 2.2 for Mbuigc '^^ 2.0 10^° M© based on a model of the COBE 
data (Zhao, Spergel & Rich 1995), r_6 ~ 1.3 for Mbuigc - 1-6 10^° Mq based on 
a revised COBE model (Han & Gould 1995), r_6 ^ 1-5 for Mbuigc =^ 1-5 10^° M© 



based on a model for the distribution of clump giants (Stanek et al. 1997) 
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r_6 '=^ 1.1 for Mbuigc ^ 1.71O^°M0 based on AGB stars (Nikolaev k Weinberg 



1997), and r_6 — 0.9 — 1.3 for Mbuige — 8 10^ M0 based on the deprojection 
of the COBE L-band data described in §2 (Bissantz et al. 1997). The latter 
paper calibrates the mass of the disk and bulge by matching the predicted gas 
velocities to the Galactic terminal velocity curve, whereas most others normalize 
their models to the velocity dispersion in Baade's window (less secure because 
of foreground disk contamination). Part of the differences between the model 
predictions is due to differences in the bulge mass. 

Thus even in bar models it has proven difficult to obtain the measured 
optical depths. One interesting result was noted by Fux (1997): in his N-body 
simulations a strong spiral arm in front of the Galactic center could add as 
much as r_6 = 0.5. The resolution of the large optical depth problem is clearly 
important. 



4. Stellar dynamics of the Galactic bulge bar 

Do measurements of stellar kinematics in low-extinction bulge fields support the 
triaxial nature of the bulge? How can they help to elucidate its structure and to 
constrain the pattern speed or aspect angle? To answer these questions requires 
dynamical models which link the stellar density distribution and gravitational 
potential to the kinematics of the stars observed. 

Before discussing existing dynamical models, it is worthwhile to note that 
Galactic kinematic observations in principle contain a very large amount of in- 
formation, (i) We can measure radial velocities as well as proper motions, (ii) 
With enough stars, velocity histograms can be obtained rather than just first 
and second moments, (iii) Distance information is available through the known 
luminosity distribution of the tracers. For the clump giant stars already dis- 
cussed, the absolute magnitude width corresponds to a distance uncertainty of 
~ 15%. The drawback of Inner Galaxy kinematic observations is, as always, 
the heavy extinction. However, if the samples are well-understood, and if it is 
possible only in a number of bulge fields to obtain an approximate measurement 
of the distribution of stars as a function of three velocities and LOS distance, 



such data will be highly constraining for the models, and thus worthwhile to 
get! 

An influential dynamical model has been that of Kent (1992). He con- 
structed an oblate-isotropic rotator model which successfully reproduced essen- 
tially all the then available velocity dispersion data in the bulge. This in itself 
is interesting because it suggests that the signatures of triaxiality in these bulge 
fields are not strong. Proper motion measurements in Baade's window (BW, 
Spaenhauer, Jones & Whitford 1992) confirmed that the velocity ellipsoid in 
BW is nearly isotropic. 

One clear signature of triaxiality is rotation on the bulge minor axis. Sofar 
the observational situation in the bulge is suggestive but not yet clear (Tyson &: 
Rich 1991, Blum et al. 1994, Izumiura et al. 1995). Another signature is a tilt 
of the velocity ellipsoid out of the (i?, z)-plane. Zhao, Spergel &: Rich (1994) 
analyzed a subsample of the proper motions stars in BW with measured radial 
velocities and found preliminary evidence for such a vertex deviation, perhaps 
evidence for the bar. 

New triaxial models for the Galactic bulge are broadly consistent with a 
range of stellar kinematic data in the bulge region. Zhao (1996) constructed a 
dynamical model for one of the Dwek et al. (1995) COBE bar models augmented 
by a central cusp. Fux (1997) analyzed N-body models whose surface density 
is consistent with the COBE K-band data after dereddening with a foreground 
screen model. Finally, Hafner et al. (1999) constructed a dynamical model by 
an extended Schwarzschild method that reproduces the deprojected luminosity 
density of BGS97 and simultaneously fits a number of kinematic constraints. The 
latter work emphasises the importance of understanding the 'selection function', 
which describes the relative fraction of stars contributing to a given kinematic 
sample, as a function of distance along the LOS. By a careful choice of sample, 
kinematic information can be obtained for a fairly narrow range of distances 
along the LOS. Such data will be the most helpful for discriminating between 
models. 

While the presently available stellar kinematics are thus compatible with 
the COBE bar models, finding unambiguous evidence for triaxiality or even 
constraining the bar parameters will require more and more specific data. 

5. Gas dynamics in the Milky Way 

It has long been suggested that the non-circular and forbidden velocities in the 
HI and CO (/, v) diagrams are the signature of a rotating bar in the inner Galaxy 
(see §1 for references). Considerable understanding can be reached with simple 
orbit-based models, because away from resonances the gas flow is expected to 
approximately follow closed ballistic orbits (e.g., Athanassoula 1992). Within 
such a model, we suggested (Binney et al. 1991) that (i) the peak at / ~ 2° and 
rapid fall-off towards larger longitudes of the Galactic terminal velocity curve 
can be explained in terms of the properties of xi-orbits in a barred bulge viewed 
nearly end-on, (ii) the dense Galactic center clouds seen in CS are moving on 
X2-orbits further in, at i? ~ 100 pc, (iii) the transition between these two flow 
regimes occurs through a cusped orbit (dust lane) shock whose trace is the 
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Figure 3. Radii of major resonances in the standard COBE (p = 20° bar, 
disk, and nucleus model of Englmaier & Gerhard (1999), for pattern speeds 
corresponding to a range of corotation radii Re = 3.0 — 4.5 kpc. 



parallelogram seen in the molecular gas observations of Bally et al. (1988), and 
(iv) the molecular ring at i? ~ 3.5 kpc is associated with the bar's OLR. 

Closed orbit models cannot take into account the effects of the strong shocks 
that form in barred galaxy gas flows (e.g., Athanassoula 1992, Englmaier & Ger- 
hard 1997). Where these shocks occur, hydrodynamic forces become important 
and so the flow may deviate strongly from ballistic orbits. Also, the Binney 
et al. analysis was based on a simple mass model for the inner Galaxy; now that 
quantitative models can be obtained from the NIR emission of old stars, the gas 
flow in the inner Galaxy can be investigated with much more confidence. Several 
recent studies have thus made a fresh attack on understanding the Galactic (/, v) 
diagrams, trying to model the non-circular motions in the central ~ 10 degrees, 
the dust lane shocks, the nature of the 3-kpc-arm, the molecular ring, and the 
Galactic spiral arms. These studies have followed complementary approaches: 
Englmaier & Gerhard (1999, EG99) have modelled SPH gas flows in the depro- 
jected COBE bulge and disk, Fux (1999) has constructed ab initio N-body/SPH 
self-consistent barred galaxy models adapted to the Milky Way, and Weiner &: 
Sellwood (1999) have studied 2D fiuids in a family of barred potentials. 

A primary question concerns the pattern speed Op, or corotation radius 
RcR, of the Galactic bar. We can estimate this as follows: (i) The 3-kpc-arm 
is associated with non-circular velocities Vnc — 50kms~^ at / = 0. From the 
results of EG99, only arms inside corotation and driven by the bar have such 
large Vnc- (ii) If the rotating potential is approximately time-independent, the 
inner edge of the molecular ring should be a lower limit to Rqr, both if this is 
indeed a ring near the outer Lindblad resonance (OLR; Schwarz 1981), and also 



if in reality it consists of tightly wound spiral arms (Dame 1993, EG99). The 
inner edge of the molecular ring can be estimated from the distribution of hot 
cloud cores (^ 4kpc, Solomon et al. 1985) and from the model of Kent et al. 
(1991) for the IRT 2.4 ^u ring (— > 3.7 kpc). Combining both arguments we obtain 
i?CR-3.5±0.5kpc. 

An independent argument from the length of the NIR bar agrees with this 
estimate. In N-body models (Sellwood & Wilkinson 1993), and from both direct 
pattern speed measurements (Merrifield & Kuijken 1995) and gas-dynamical 
modeUing (Athanassoula 1992), it is found that i?cR — 1-0 — 1.2 x i?bar- With 
i?bar ^ 3.2 ± 0.3 kpc (§2, Preudenreich 1998) this results in i^cR - 3.5 ± 0.4 kpc 
(not including possible systematic errors in locating the end of the bar). 

However, in self-consistent N-body/SPH models Fux (1999) has observed 
time-dependent gas flows that never approach a quasi-stationary configuration. 
The time-dependence appears to be driven partly by a central sloshing of the 
bulge with amplitude of several 100 pc, which is not well understood, and partly 
by the interaction of the bar with the outer stellar and gaseous spiral arms. 
If the Milky Way is a similarly time-dependent system (see Morris & Serabyn 
1996), argument (ii) above does not necessarily hold because then the spiral 
arms in the molecular ring may penetrate into corotation. Fitting his models to 
the observed (/,f) diagrams, Fux (1999) finds a larger -RcR = 4.0 ± 0.5 kpc. 

Fig. 3 shows the radii of the major resonances for the combined range 
of corotation radii i?cR = 3.0 — 4.5 kpc and a bulge, disk and nucleus mass 
distribution based on the 20° model of BGS97 with constant L-band mass-to- 
light ratio. The absolute scaling of the ordinate is approximate, depending on 
the detailed model fit to the terminal velocity curve. The position of the OLR 
is predicted to be at 5.3 — 6.7 kpc, so is well inside the solar radius even for the 
slowest pattern speed. The circular velocity at these radii from the luminous 
matter only is still close to 220kms~^, so adding a dark halo near and beyond 
the solar radius does not change the position of the OLR significantly. Thus the 
influence of the bar on the dynamics of the disk near the Sun should be small. 

Is this confirmed by more detailed simulations? Fig. 4 compares model gas 
flows in the full = 20° COBE potential and in the same model modifled by 
axisymmetrizing the disk for radii R > 3.5 kpc. The two gas flows are dramat- 
ically different; the bar alone is unable to drive signiflcant spiral arm structure 
in the gas beyond R ~ 4 kpc - the structures seen in the right panel of Fig. 4 
are merely gas concentrations but no real shocks; this is not a resolution effect. 
Also, the two arms beginning above and below the bar in the left panel of Fig. 4 
are absent in the right panel. The bar-only model would be in conflict with 
the tangent point observations indicating the presence of strong spiral arms at 
I = —50°, —30°, 30°, 50°, and the abundance of warm CO clouds found near 
/ = 25 — 30° (Solomon, Sanders & Rivolo 1985) indicating a spiral arm shock 
in this region. In the full COBE model, the spiral arms beyond R ~ 4 kpc 
are driven by the mass corresponding to the NIR luminosity concentrations at 
^ 3 kpc down the bar minor axis (see §2). The conclusion is therefore that these 
concentrations cannot just be caused by luminous supergiant stars, but must 
contain significant mass. Most likely, they are due to incorrectly (because of the 
employed symmetry assumption) deprojected spiral arms. If this interpretation 
is correct, the dynamics between the bar's corotation radius and the Sun is dom- 
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Figure 4. Left: Quasi-stationary gas flow in the eightfold-symmetric 
COBE bulge and disk potential, deprojected for </) = 20°, with an added 
nuclear cusp. Corotation is at 3.4 kpc, and the Sun is to the left along the 
dashed line. Right: Gas flow in the same model, but at lower resolution 
and with the density multipoles with ?Ti 7^ set to zero outside 3 kpc. From 
Englmaier & Gerhard (1999). 



inated by the Galaxy's spiral arms. These might well have a pattern speed of 
their own (Sellwood h, Wilkinson 1993), and then they and not the bar would 
dominate resonances near the Sun. 

Outside -RcR the induced non-circular motions are fairly small. Fig. 5 
shows predicted terminal velocity curves from several gas flow models in the 
gravitational potential of the COBE bar and disk, to which in some cases a 
dark halo potential has been added. In these models the bar rotates with a 
pattern speed corresponding to -RcR — 3.4 kpc (see EG99). Overall, the models 
match the observed HI and CO terminal velocities rather well. The two main 
discrepancies are in the central peak, where there are both resolution problems 
and uncertainties in the potential, and in a region around I ~ —20°, where the 
disk potential from the deprojected NIR data may be seriously in error (see 
§2). Out to -R > 5 kpc, where the dark halo begins to contribute, the observed 
terminal velocity curve is consistent with the radial distribution of old stellar 
mass inferred from the COBE data. 

These gas flow models give rise to (/, v) diagrams which qualitatively re- 
produce a number of features seen in the observed (/,f) diagrams (see EG99 
and the paper by Englmaier & Gerhard in these proceedings). The four-armed 
spiral structure outside corotation (see Fig. 4) is quantitatively consistent with 
the observed directions to the five main spiral arm tangents at |/| < 60°. The 
3-kpc-arm is identified with one of the model arms emanating from the ends of 
the bar and extending into the corotation region (below the bar in Fig. 4). The 
cusped-orbit shock transition channels the gas onto an inner disk on a;2-orbits 
which in the simulations of EG99 has radius ~ 150 pc and rotation velocities 
~ lOOkms"^, approximately as inferred from CS observations (Bally tt al. 1988, 
Binney 1994). The models fail, however, in accounting for the magnitude and 
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Figure 5. Terminal velocity curves of several gas flow models scaled to the 
observed HI and CO terminal velocities (symbols; see EG99 for references). 
Models are similar to that in Fig. 4: without halo (dot-dashed); including a 
dark halo with asymptotic circular speed Vc = 208kms~^, highest resolution 
model (full line); including a halo with Vc = 220 km s^^ (dashed line). 



extent in longitude of the forbidden velocities seen in low-intensity HI observa- 
tions. Weiner & Sellwood (1999) have specifically modelled these HI data, and 
find that this problem can be resolved if the Galactic bar is inclined by a larger 
angle, (j) ~ 30°, as seen from the Sun. The models of Fux (1999) show, however, 
that extended forbidden velocity regions can occur also for smaller bar angles if 
the model is asymmetric and no longer quasi-stationary. Thus there may not 
be a conflict with the smaller cp favoured by the clump giant distribution and 
the large microlensing optical depth. 

The non-stationary gas flow models in the N-body and SPH simulations of 
Fux (1999; see also his paper in these proceedings) give the best match sofar to 
several features in the observed {I, v) diagrams: the 3-kpc arm, the connecting 
arm, the 135kms~^-arm, and the molecular ring. They fail in that there is no 
spiral arm corresponding to the / = 50° tangent, and they assert that the inner 
disk on a;2-orbits rotating at ~ 200kms~^ extends to i? ~ Ikpc, distinct from 
the dense clouds seen in CS which extend to i? ~ 2° and appear to rotate with 
~ lOOkms"^ on presumably j;2-orbits (Binney 1994). The issue of a large X2- 
disk is also relevant for interpreting the kinematics of the bar-driven inner arms. 
The interpretation of the 3-kpc arm as one of the lateral arms emanating from 
the end of the bar is similar in the models of EG99 and Fux (1999), although 
the morphology of the transition from this arm into the dust lane shock is 
somewhat different, probably due to differences in the gravitational potential 
and differences between quasi-stationary and time-dependent flows. 
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Sevenster (1997) argues that in reality the 3-kpc arm is part of an inner 
ring surrounding the bar, such as observed in some other galaxies. Her strongest 
observational result supporting this interpretation is (in my view) the alignment 
of a small group of 100-300 Myr old OH/IR stars coincident with the 3-kpc arm 
in the HI {l,v) diagram. The rotation period at 3 kpc is ^ 100 Myr, so these 
stars appear to have remained near their birth sites for 1-3 rotation periods. 
Much depends on the precise ages of these stars, and on the pattern speed of 
the 3-kpc arm. From the above this will be similar to the bar's pattern speed. 
Stars at 3 kpc radius then leave the arm on a time-scale Prcl = 27r/[r2(3kpc) — 
ftp] ~ 27r/0(3kpc) X [1 — 3kpc/i?cR]^^ (for flat rotation curve). With Rqr = 
[4.5 kpc, 3.5 kpc] Pj-ei is [3x,7x] the orbital period, and for Ren = 3kpc these 
stars would not leave the pattern at all. 

The now existing models for the Galactic gas flow certainly do not yet 
provide a quantitative fit to all the major features in the observed {l,v) dia- 
grams. However, with their help many of these features are beginning to make 
sense in the framework of a barred Milky Way. It is encouraging that two ap- 
proaches, one starting from the COBE observations, the other from the evolution 
of self-consistent N-body systems, seem to converge to a similar picture when 
constrained by the HI and and CO {l,v) diagrams: a corotation radius about 
half way between the Sun and the Galactic Center, a strong, probably 4-armed 
spiral pattern between i?cR ^-nd the solar radius, and non-circular gas motions 
inside RcR which take the gas along the inner arms near the ends of the bar 
through the dust-lane shocks to the nuclear disk on X2-orbits. 



6. The mass of the bulge and disk, the Oort hmit, and the 
distribution of dark matter in the Galaxy 

The observed terminal velocity curve (Fig. 5) is consistent with a maximum NIR 
disk (and bulge) model out to |/| ~ 45° or i? ~ 5.5 kpc ~ 2i?/), for an exponential 
scale length of the disk of R^ ~ 2.5 kpc (see §2). This maximum disk model 
predicts a surface mass density of S© = 44 — 49M0/pc^ near the Sun at R0 = 
8 kpc for a local circular velocity of Vc,q = 208 — 220kms~^. The local surface 
density of 'identified matter' is 48 ± 8Mq/pc^, while a combined constraint 
from the total surface density within z = ±1.1 kpc and the rotation curve is 
48 lb 9Mq/pc^ (Kuijken & Gilmore 1991). Both are consistent also with the 
K-giant analysis of Flynn & Fuchs (1994). That these numbers approximately 
agree lends support to the conclusion that the Galaxy indeed has a near-maximal 
disk. Similarly, the high microlensing optical depth to the bulge argues for a 
high bulge and disk mass, but a quantitative argument is presently difficult to 
make (see §3). Compared to earlier analyses, the main difference is the short 
disk scale-length (see also Sackett 1997) and the extra light around 3 kpc (see 
BGS97) - the Sun is well beyond the maximum in the rotation curve from the 
NIR luminous matter only. 

The mass of the NIR disk inferred from the terminal velocities must include 
some but may not include all of the mass in the Milky Way's gas disk (local sur- 
face density ~ IOMq/pc^; see Dame 1993) and thick disk (local surface density 
~ 9M0/pc^ and scale length i?th — 4.5 kpc. Beers & Sommer-Larsen 1995, 
Ojha et al. 1996), because the radial density distributions of these components 
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are different from the old NIR disk. For the following I have assumed that the 
total visible mass in the Galaxy is given by the mass of the COBE NIR bulge 
and disk and the thick disk only, assuming that the gas disk is already accounted 
for by the fitted mass in the NIR disk. Since the HI distribution has a larger 
scale-length than the stars, this approximation might slightly overestimate the 
final halo contribution. 

Computing the gravitational forces from these two components, and fitting 
a cored spherical halo so as to make the Galaxy's rotation curve approximately 
flat at R> 6kpc with Vc = 220kms~^, results in a halo core radius Re ^ 15kpc. 
Integrating the surface density of this halo between z = ±1.1 kpc gives S/i,i.i = 
17 Mq/ pc^. Adding this to the surface density of the old NIR disk and the thick 
disk, the total is Snir + ^th + S/i^i.i = 70 — 75M0/pc^, whereas the measured 
total Si.i = 71 ± GMo/pc^ (Kuijken & Gilmore 1991). 

This very good agreement is clearly better than one expects, given the 
various uncertainties in both numbers, but it shows that a maximum NIR disk 
model and a near-spherical halo with a large core radius provide a natural 
explanation for both the observed terminal velocity curve and the measured 
surface density near the Sun. Note that for smaller values of the Galaxy's 
asymptotic rotation velocity, the required amount of halo would be reduced; for 
Vc = 180 km s^"*^ the terminal velocity curve can be fit within the errors without 
any added halo (at the cost of a falling rotation curve) . A more detailed analysis 
is clearly worthwhile. 

7. Concluding remarks 

The work reviewed here confirms that the Galaxy contains a central bar, seen 
most clearly in the NIR, the source counts and the gas kinematic observations. 
The corotation radius of the bar is about half-way between the Sun and the 
Galactic Center. There is still some uncertainty about the orientation (most 
likely cf) = 20 — 25°, but possibly = 15 — 35°) and axial ratios (roughly 3:2:1) 
of the bar, and about the relative importance of the disk in the central few kpc. 
The short disk scale-length in the NIR is important; combined with the Galactic 
terminal velocity curve and the Oort limit it implies a near-maximal disk and 
a fairly large halo core radius. 

Is there a Galactic bulge besides the bar? This has not been conclusively 
answered. While the clump giant stars appear to belong to a strongly barred 
component, the NIR-luminosity could be less strongly barred, and the RR Lyrae 
stars, presumably part of the stellar halo, appear to be unbarred (Alcock et al. 
1998). In external galaxies, kpc-scale central bars and nuclear bulges may co- 
incide (MollenhofF, private communication; see also Seigar & James 1998). In 
M94, the vertical velocity dispersion in the nuclear region is higher than that 
in the bar (Mollenhoff et al. 1995); this suggests that in the Milky Way the 
question might be answered by studying proper motions as a function of age 
and metallicity in low latitude bulge fields. 
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